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We implement an effective operator formalism for general one- and two-body operators, obtaining 
results consistent with the no-core shell model (NCSM) wave functions. The Argonne V8' nucleon- 
nucleon potential was used in order to obtain realistic wave functions for ''He, ^Li and '^^C. In the 
NCSM formalism, we compute electromagnetic properties using the two-body cluster approximation 
for the effective operators and obtain results which are sensitive to the range of the bare operator. 
To illuminate the dependence on the range, we employ a Gaussian two-body operator of variable 
range, finding weak renormalization of long range operators (e.g., quadrupole) in a fixed model 
space. This is understood in terms of the two-body cluster approximation which accounts mainly 
for short-range correlations. Consequently, short range operators, such as the relative kinetic energy, 
will be well renormalized in the two-body cluster approximation. 

PACS numbers: 21.60.Cs, 23.20.-g, 23.20.Js 



I. INTRODUCTION 

The no-core shell model (NCSM) is a successful ab ini- 
tio method used to compute properties of light nuclei, 
starting from realistic two- and three- nucleon interac- 
tions. In this approach, all nucleons are active, but the 
Hilbert space available is a finite model space. In or- 
der to take into account the effect of the configurations 
outside the model space, we employ effective interactions 
obtained by means of an unitary transformation 
aimed to reproduce the low-lying spectra in our finite 
model space. Large basis NCSM calculations have been 
successful in the description of energy spectra in light 
nuclei 0, 0, 0]. However, the E2 transition strengths, 
powerful tests of the theoretical wave functions, are usu- 
ally underestimated 0,01 • To date, they were obtained 
with the bare operators, and one expects that a renor- 
malization in the same fashion as for the interaction will 
give a better description. The only such effective opera- 
tors computed so far have been the effective point-proton 
radius 0], relative kinetic energy 0, and the nuclcon- 
nucleon (NN) pair density 0, U • An earlier application 
to electromagnetic operators was in a restricted model, 
where the effective quadrupole charges have been shown 
to have the values expected in a single harmonic oscilla- 
tor shell, while the effective two-body contributions were 
found to be small Q. 

The purpose of this paper is to present the formalism 
for general one- and two-body effective operators consis- 
tent with the effective interactions of the ab initio NCSM. 
Section HTl details this approach, particularly with regard 
to the treatment of states outside the model space. In 
Sec. Illll we apply the formalism to different one- and two- 
body operators. The renormalization is found to depend 



upon the range of the operator, weak for the E2 opera- 
tor, and strong for the relative kinetic energy. Because 
the two-body cluster approximation used in this work 
integrates out mainly the short-range part of the inter- 
action, short-range operators, such as the relative kinetic 
energy, arc strongly renormalized, while long-range op- 
erators, e.g., quadrupole, are weakly renormalized. By 
using a Gaussian two-body operator of variable range, 
we show how much the renormalization depends on the 
range of the operator. We present our conclusions in Sec. 



II. THEORETICAL OVERVIEW 

We start with a system of A particles, interacting 
through the intrinsic Hamiltonian 
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where m is the nucleon mass, and V^^ the bare NN in- 
teraction, such as the Argonne potentials in coordinate 
space hfl or the non-local CD-Bonn pd] |. While realis- 
tic three-body forces have been shown to be important 
in obtaining the nuclear spectra [T^ IT^ and in de- 
scribing electromagnetic and weak form factors [T5| . we 
consider only two-body interactions for simplicity. 

The intrinsic properties of the system are not affected 
by the addition of the CM harmonic oscillator (HO) 
Hamiltonian, but by casting the new Hamiltonian in the 
form 



2 



H'l = Ha 



p2 



1 



2mA 2 



E 

i=l 
A 

E 



— + -. 

2m 2 



mf}2,f 



'J' 2A 



A A 

= E + E ' 

z— 1 ?>j^l 



(2) 



and applying the unitary transformation on the new 
Hamiltonian, we improve the rate of convergence of the 
solution of the A-body problem in smaller model spaces. 
As we subtract the CM term in the final many-body cal- 
culation, it docs not introduce any net influence on the 
converged intrinsic properties of the many-body calcu- 
lation. Furthermore, this addition and subtraction does 
not affect our exact treatment of the CM motion. This 
procedure introduces a pseudo-dependence upon the HO 
frequency and the two-body cluster approximation de- 
scribed below will sense this dependence. In the largest 
model spaces, however, important observables manifest 
a considerable independence of the frequency and the 
model space limit. 

A transformation, which accommodates the short 
range correlations by means of an antihermitian oper- 
ator 0,1113, produces an effective Hamiltonian H. given 

by 
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Note that even if the original Hamiltonian contained just 
one- and two-body terms, the operator S and the trans- 
formed Hamiltonian Ti, contain up to A-body terms. Ob- 
taining the exact operator S is equivalent to solving the 
initial problem, which would make the transformation 
impractical. Wc will discuss below, in some detail, ap- 
proximations which allow us to solve the A-body system; 
for now, the derivation is exact. 

The purpose of transformation is to preserve the 
solutions of the original Hamiltonian when one reduces 
the infinite dimension of the Hilbert space to a finite 
value which is numerically tractable. We achieve this by 
splitting the full space associated with the A-body system 
into P, or model space, and Q, the excluded space; the 
decoupling condition 



QHP = 0, 
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in addition to the requirements PSP = QSQ = 00, en- 
sure that correlations left out by restriction to the model 
space P are properly taken into account for a subset of 
the exact eigenstates. Formally, the operator S can be 
written Q by means of another operator w as 



5* = arctanh(ti> — oj^). 
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where the new operator fulfills QluP = lo. The energy- 
independent effective Hamiltonian in the model space P 
becomes 



Fe// = PUP = 



P + Puo^Q ^,^ P + QuP 
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and, analogously, any arbitrary observable can be trans- 
formed to the P space as 
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In order to compute effective operators one needs to 
know uj. This operator connects eigenvectors in P to 
vectors in the Q space. A simple way to compute w is Q 
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with \ap) and jag) the basis states of the P and Q spaces, 
respectively; |fc) denotes states from a selected set K. of 
eigenvectors of the Hamiltonian in the full space \ k) = 
Ek\k), and {ap\k) is the matrix element of the inverse 

overlap matrix {ap\k), that is J2api^'\'^p)i'^p\^) ~ ^kk'- 
Note that the dimension of the subspace /C is equal with 
the dimension of the model space P. 

As noted before and seen explicitly in Eq. (jHJ, in or- 
der to solve for ui one needs the solution of the A-body 
problem, i.e., the eigenvectors \k), which is the final goal. 
Therefore, we introduce the cluster approximation. This 
consists in finding uj for the a-body problem, a < A, and 
then using the effective interaction thus obtained for solv- 
ing the A-body system. This approximation introduces 
a real dependence of the oscillator parameter il, and the 
solution to this problem is to search for a range of val- 
ues over which the results are weakly n dependent. There 
arc two limiting cases of the cluster approximation: first, 
when a — > A, the solution becomes exact; a higher-order 
cluster is a better approximation and was shown to in- 
crease the rate of convergence Second, when P — > 1, 
the effective interaction approaches the bare interaction; 
as a result, the cluster approximation effects can be min- 
imized by increasing as much as possible the size of the 
model-space size. 

In this work, we present results obtained at the two- 
body cluster level. Under this approximation, the trans- 
formation writes as 



5* w ^ Sij , 
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with Sij = arctanh(cjy — i^lj)- Applying the operator 
identity 



e-^Oe" =0 + [0,S] + - [[0,S],S' 



(10) 
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to transform a general one-body operator O^^^ = 
SiLi 0116 obtains 

A A 

= o(i)+ ^ [0,+0„S,j]+Y}\0^+0„S,,],S,,]+..., 

i>j=l i>j 

(11) 

where we have retained only the one- and two-body 
terms, neglecting higher body contributions, such as 
[0^,5^7;], with i ^ j and i ^ k. Resummation of the 
commutators yields 

POeffP = PY.O^P 

i 

A 

+P [e"^'^ {0^ + Oj) e^- - (O, + O,)] P.{12) 

i>j=i 

Analogously, for a general two-body operator 

A 

POeffP = P J2 e-^*^Oye^-P, (13) 
i>j=i 

and, in particular, the effective Hamiltonian derived from 
Eq. (O is given by 

A 

PHeffP = Pj2h,P 
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We emphasize that in the two-body cluster approxima- 
tion the explicit decoupling condition in Eq. Q) is ful- 
filled for the two-body problem now: 

Q2H^'^P2 = 02e-^^^(/ii + h2 + vi2)e^'^P2 = 0, 

where P2, Q2 refer to the corresponding projection op- 
erators for the two-particle system. Condition is, in 
general, violated for the A-body problem, but the errors 
become smaller with increasing the model space. The 
dependence upon A due to the addition of the CM term 
in Eq. jSJ is kept in vi2. 



Finally, note that even if initially one starts with an 
one-body operator, in the two-body cluster approxima- 
tion, the effective operator will generally have irreducible 
two-body matrix elements. 



III. RESULTS AND DISCUSSION 

For details regarding the procedure to obtain the ef- 
fective interaction for a system of A nucleons in the 
two-body cluster approximation, we refer the interested 
reader to previous work, e.g., Ref. 0. We note that in 
this paper the effective interaction and the transforma- 
tion LO are obtained in the relative system of two particles, 
in a large HO basis. The Q space is chosen to be a few 
hundred fi^ excitations in order to obtain an exact solu- 
tion to the two-body Schrodinger equation. Due to the 
rotational symmetry, we formulate the problem in two- 
nucleon channels with good total spin s, total angular 
momentum j, and isospin t, reducing drastically the di- 
mensions involved, when performing the summation over 
the states in the Q space in Eq. ^ The same proce- 
dure can be applied to operators which can be analyti- 
cally expressed in terms of relative and CM coordinates 
of pairs. Therefore, we develop a convergence procedure 
which works for general one- and two-body operators. 

As shown in the previous section, the corrections to 
general operators arc given by Eq. (|12|) for one-body op- 
erators, and Eq. for two-body operators. Suppose 
the operators are given in the single-particle representa- 
tion, so that in order to compute contributions of the 
form exp(— S'ij)(Oi H- Oj ) exp(S'ij) for one-body opera- 
tors [and correspondingly exp(— S'ij)Oy exp(S'y) for two- 
body operators] by means of Eq. Q , one needs to trans- 
form either the operator to the relative system or the 
transformation lo to single-particle representation. The 
two procedures, however, give the same result, and the 
numerical burden is likely comparable. 

Going back to Eq. |(7J), note that in the computation 
of effective operators, the number of two-body matrix el- 
ements involved in the summation over the Q states be- 
comes numerically intractable in the single-particle repre- 
sentation. For example, in terms of the matrix elements 
in the relative coordinates, the transformation uj in the 
single-particle picture is given by 
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where \'niliji,n2l2j2] Jt) refer to two-body states in the 
Q2 space, and 1713/3^3, 714^4^4; Jt) refer to two-body states 
in the P2 space. We have employed the same notations 
as in Ref. [13I for the Talmi-Moshinsky transformation. 

Analyzing Eq. (|15|l . we note that the quantum num- 
bers for the CM states are restricted by the model space, 
2N + L < Np, with Np fixed by the size of the model 
space. The relative {n, I) states are restricted by the 
model space so that 2n + I < Np. Using the en- 
ergy conservation in the Brody-Moshinsky brackets, one 
also obtains a restriction for the single-particles states, 
2??3 + I3 + 2n4 -I- ^4 < Np. However, the states (n' , I') run 
over the excluded space and, since 27t.i -I- Zi 4- 2n2 + I2 is 
not restricted by the model space, the number of possible 
pairs becomes numerically intractable. Hence, we restrict 
the relative states in the excluded space by the condition 
Np < 2n' + 1' < Nq , and observe convergence by increas- 
ing Nq . (To give the reader an idea about the dimensions 
involved, we note that for Nq = 28 for a p-shell nucleus, 
the number of lu matrix elements, taking into account 
the possible symmetries, in a 2hQ model space defined 
by Np = 4, is 413,163; subsequently, the number of re- 
duced matrix elements involved in the transformation of 
a tensor operator with J = 1, T = 1 is 7,216,180.) This 
procedure was successfully tested for the deuteron in a 
restricted space, where we have shown that, for general 
tensor operators, the matrix elements obtained with ef- 
fective operators approach the values obtained with bare 
operators in the full space [T^ . 

In order to test the convergence procedure in a realis- 
tic model space, one could consider two-body operators 
which depend only upon the relative coordinates, which 
can be renormalized similarly to the Hamiltonian. One 
such operator is the relative kinetic energy, whose effec- 
tive matrix elements can be computed by simply replac- 
ing in the effective interaction code the bare Hamilto- 
nian with the bare kinetic energy after the transforma- 
tion LU has been determined. The results summarized in 
Fig. n show that the expectation value computed with 
the approximate effective operator converges toward the 
effective value computed including all the states in the 
excluded space. While the convergence rate might look 
slow, we would like to point out that even for small values 
of Nq we obtain reasonable renormalization comparing 




ii OOShQ - 



bare 12 18 24 30 Nq^co 
Nq 

FIG. 1: (Color online) The expectation value of the relative 
kinetic energy on the ground state using realistic wave func- 
tions for *He in different model spaces. We show how the 
ground state expectation value changes when one increases 
the number of states in the Q-space (circles). For comparison, 
we include the results obtained with bare operator (squares) 
and effective operator including all states in the excluded 
space (diamonds). 



to the full space renormalization, in contrast to the in- 
teraction, where several hundred hfl excitations are usu- 
ally necessary. Note that larger model spaces require less 
renormalization with respect to the bare operator, as ex- 
pected. Also, the effective expectation values are similar 
to previous results , although the latter have been ob- 
tained with wave functions computed in the three-body 
cluster approximation in IGhfl. This can be explained by 
the character of the kinetic energy operator which is zero 
range, so that it is well renormalized at the two-body 
cluster level. 

Our first applications of the formalism to one-body op- 
erators are to quadrupole and Ml transition operators. 
We present in Fig. |21 the 2h,^l model space results for 
B{E2;2+0 ^ Oj^O) in ^^C and B{Ml;l+0 0+0) in 
^Li. The shell model calculations have been performed 
using the many-fermion dynamics code [T^ . Following 
the procedure described above, we compute the effective 
operator by means of Eq. ((T^ by adding from the Q 
space two-body matrix elements (27?.' -|-/') < Nq, and in- 
creasing Nq . Note that for E2 transitions we expect the 
biggest contribution to come from ATiil, as the E2 opera- 
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FIG. 2: (Color online) 2TiQ model space results for B{E2) 
in ^^C and B(M1) in ^Li using effective interaction derived 
from AV8' potential. We show the results as a function of the 
dimension of the Q-space included. For comparison, we in- 
cluded the values obtained using the bare operators (squares) 
and the experimental values (diamonds). The dashed lines 
represent an interpolation of the strengths obtained with ef- 
fective operators. 



tor connects across two shells. Figure |2 shows, however, 
that B{E2; 2+0 ^ O^^O) remains essentially flat, at the 
same value as the one obtained with the bare operator, 
and at about half the experimental strength. We have 
included the experimental B{E2) for reference. What 
we do not have is the exact results for the full Hilbert 
space which would provide the ultimate comparison with 
experiment. Instead, we present the convergence with in- 
creasing Nq at the 2-body cluster level as shown in Fig. 
121 Our achievement here is to include the contribution 
for renormalization at the two-body cluster level which 
turns out to be rather small. The large residual discrep- 
ancy between our extrapolated result and experiment is 
due to the combination of residual effective three-body 
effects and to the neglect of genuine three nucleon poten- 
tials. However, we argue below that the result obtained 
with only two-body interactions should be significantly 
closer to the experiment. 

For the Ml, we observe a small correction in the transi- 
tion strength, albeit the correction is slightly larger than 
the discrepancy between theory and experiment. For 
other Ml transitions, however, we observe, in general, 
the same minimal effect as for quadrupole transitions. 
Such a small correction can be more easily understood, 
as Ml does not connect the model space with the com- 
plementary one; increasing the model space introduces 
enough correlations in the wave function so that using 
the bare operator gives B{M1) values close to experi- 
ment. 



Despite the successful test in the case of the relative 
kinetic energy, one can still ask the question whether the 
convergence procedure is faulty and if including all the 
states in the excluded space in calculations would not 
improve the E2 transition strengths. Because, in gen- 
eral, the electric multipoles can be written as sum of 
terms which factorize the CM and relative contributions 
for each pair of nucleons, we can include all the states in 
the Q-space. (The corresponding transition strengths are 
analogous to the Nq — > oo values in Fig. ^) In the Ap- 
pendix, we identify the CM and relative contributions for 
each pair to the E2 operator. In this context, we further 
show that the difference between the usual E2 operator 
used in shell model calculations, and the translationally 
invariant one is a term which involves only contributions 
from the CM of the A-body system. The latter is ir- 
relevant when employing translationally invariant wave 
functions or many-body wave functions which factorizes 
exactly in intrinsic and a Ohfl CM contributions, as in 
our case. 

In TableHl we present the B{E2) values for quadrupole 
transitions in a 2hQ model space in ^Li, and compare 
them with the experimental strengths, where available. 
In these calculations we include in the summation all the 
states in the Q space. Because, as noted before, we limit 
ourselves to small model spaces, the transition strength 
obtained with bare operators are far from the experi- 
mental values. As noted before, part of the discrepancies 
can appear because of the missing three-body forces in 
our calculations. However, we argue that the theoreti- 
cal strength obtained in the model which includes only 
two-body forces should be closer to experiment. Thus, in- 
creasing the model space up to lOhfl reduces significantly 
the discrepancy between theory and experiment for the 
1+0 3+0 transition in Tabled A similar trend is ob- 
served in larger model spaces accessible today, which sug- 
gests that indeed the strength in the full space should be 
significantly closer to experiment even in the absence of 
three-body interactions. Note, however, the large differ- 
ence between the strengths obtained with bare operators 
in 2hrt and lOhft. Because of the convergence properties 
discussed for the effective interaction in the previous sec- 
tion, which can be extended to other operators, one ex- 
pects that the bare value in the larger model space would 
be closer to the exact theoretical value than the one in the 
smaller space. Therefore, one may conclude that at least 
the value in the 2hn space is far from the correct value, 
and that the effective operator should have a significant 
impact in this model space. Table shows the opposite 
result, i.e., very little difference between the transition 
strengths obtained with bare and effective operators. In 
contrast, for the deuteron, where the two-body cluster 
provides the exact solution, the bare quadrupole opera- 
tor in 4hil gives 0.179 e fm^ for the quadrupole moment, 
while the value of 0.270 e fm^, described by the AV8' 
potential, is obtained using the corresponding effective 
operator. Overall, the difference between the bare oper- 
ator results in the 2hrt and lOhil model spaces, coupled 
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FIG. 3: (Color online) Left panel: renormalization of the ground state expectation value of the relative Gaussian operator 
using realistic wave functions for ''He, as a function of the range of the operator for Ahft (circles), GTlQ, (squares), and 8TlQ 
(diamonds). Right panel: expectation values in different model spaces for Gaussian operators of selected ranges. All the states 
in the Q-space are included. 



TABLE I: B{E2) values, in fm", for ''Li computed with 
bare and effective operators in a 2hQ model space {h^l = 13 
MeV). All states in the Q space are included (see text for 
details). For comparison we have included also B{E2) values 
obtained in lOhil space with the bare quadrupole operator 
Q. Experimental values are from Ref. ^y|. 
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with the small renormalization at the two-body cluster 
level, indicate there are sizable effective multi-body in- 
teraction effects needed to correct the 2hn B{E2) value. 

Because the two-body cluster approximation accom- 
modates the short range correlations, one can expect that 
such an approach might not be well-suited to renormal- 
ize the E2 operator, which is infinite range. To illustrate 
the importance of the range of the operator for renor- 
malization at the two-body cluster level, we consider a 
Gaussian two-body operator of range ap 

0(ri, f2) = Co exp (^-^H-p)!^ , (16) 

with Cq chosen so that 

Co J dr exp ^ 2^ = 1- 

While this is not a realistic observable, one can often 
expand realistic operators as sums of Gaussians, so this 



could be used to estimate the renormalization of differ- 
ent contributions. However, in this paper, the only pur- 
pose of this example is to illustrate the dependence of 
the renormalization upon the range of the operator. We 
define the renormalization as ((Og//) — {Obare)) / {Obare) , 
and in Fig. OJa) we summarize the results using the real- 
istic ground-state wave function for "^He. At small ranges, 
the expectation value computed with the effective oper- 
ator is significantly different from the one obtained using 
the bare operator. However, when the operator becomes 
longer range, the renormalized value becomes nearly in- 
distinguishable from the bare value. But in the absence 
of an exact full space expectation value of this operator 
one cannot tell whether or not the absence of renormal- 
ization for long range operators is a consequence of the 
fact that the bare operator is already close to the exact 
result. However, in Fig. Efb) we show the dependence of 
the ground-state expectation value of the operator upon 
the size of the model space for selected short- and long- 
range operators. For the short-range operator (aq = 0.2 
fm) , the expectation value obtained with the bare opera- 
tor varies with the space, while the one obtained with the 
effective operator is flat, suggesting that indeed the cal- 
culation is converged in this case. On the other hand, the 
expectation value of the long-range operator (oq = 1 fm) 
presents about the same dependence on the model-space 
size with both bare and effective operators, therefore sug- 
gesting that the effective operator at this range is as poor 
an approximation to the exact result of the full space as 
the bare operator - a feature reminiscent of our B{E2) 
results above. That is, while we do not expect realis- 
tic operators to behave exactly as the Gaussian operator 
used here, we believe that this example offers a qualita- 
tive understanding of the very weak renormalization of 
the quadrupole operator. 
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IV. SUMMARY AND CONCLUSIONS 

We have computed the effective operators consistent 
with the ab initio NCSM effective interactions. In the 
present investigation we hmit ourselves to solutions ob- 
tained in the two-body cluster approximation. The main 
purpose of this paper was a qualitative investigation of 
the effective operators rather than a highly accurate de- 
scription of the experimental data; this is the reason we 
have considered limited shell model spaces for ''He, ^Li 
and ^^C nuclei. Such small model spaces offer an excel- 
lent testing ground for effective operators, for one expects 
that the smaller the model space, the more significant the 
renormalization. 

For general operators one cannot apply the same cal- 
culational procedure as for the Hamiltonian; for these 
we have developed a convergence procedure in which we 
add correlations from the excluded space a few shells at 
a time. This procedure can be applied to arbitrary one- 
and two-body operators. We have tested this procedure 
in the case of the relative kinetic energy, where we ob- 
tained slow convergence as the number of Q-states in- 
cluded in the renormalization is increased. 

We found that the quadrupole operator is very weakly 
renormalized at the two-body cluster level. We suggest 
that this is a consequence of the character of multipole 
operators which have infinite range. To substantiate this 
hypothesis, we have shown, using a Gaussian operator 
that the renormalized operator is quite sensitive to the 
range of the operator. Thus, the shorter the range, the 
stronger the renormalization, with only a small effect for 
long range operators. Therefore, in order to describe long 
range operators one needs to go beyond the two-body 
cluster approximation. As a general caveat, any trun- 
cation of the space could induce effective operators with 
non-negligible higher-body correlations. This result is in 
accord with previous findings of tests in restricted models 
for the double-/3 decay operator . It should be noted 
that earlier calculations for ^Li, which obtained signif- 
icant effective quadrupole charge renormalization, were 
based on large-basis NCSM calculations, which were then 
explicitly truncated into a Ohfl space and fitted to one- 
plus-two-body quadrupole operators 0. By construc- 
tion these calculations contained all correlations up to 
six-body due to the truncation and, hence, yielded the 
large effective charge renormalizations that are found ex- 
perimentally. Techniques for including these high-body 
correlations in our calculations are under investigation. 

Nevertheless, even at the two-body cluster level, there 
are cases where the renormalization is significant, as we 
have shown for the relative kinetic energy. Future work 
will investigate electromagnetic processes, where one ex- 
pects the two-body cluster renormalization to play an 
important role at large momentum transfer. 
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APPENDIX: THE TRANSLATIONALLY 
INVARIANT QUADRUPOLE OPERATOR 

In this Appendix we revisit the quadrupole operator, 
identifying the intrinsic and CM contributions. 

We start by rewriting the single-particle isoscalar 
quadrupole operator in a two-body form, that is 

E2 = J20{n)^j^ J2 (O(r0 + O(r-)), (A.l) 

i j>j = l 

where 0{ri) = ejsrfY2{fi), with Y2 the spherical har- 
monics of rank 2 and e/s the isoscalar charge. For a pair 
of particles {i, j) we define the relative = fi — fj and 
the CM Rij = {ri +rj)/2 coordinates. Using Eq. (35) in 
Ref. jl^l , one can show that 

Oin) + 0{r,) = \o{n,) + 20(4). (A.2) 

This separates for a pair of particles the relative and 
CM contributions to the quadrupole operator. The Lee- 
Suzuki transformation affects only the relative coordi- 
nates, so that in this form one can apply a similar pro- 
cedure as for the renormalization of the Hamiltonian. 
Note, however, that because the quadrupole operator is 
a tensor operator of rank 2, it can couple different chan- 
nels with different total angular momentum 7, unlike the 
Hamiltonian. Hence, the renormalization procedure for 
the E2 multipole is more involved than for the Hamilto- 
nian. 

We can show in general that any electric multipole op- 
erator can be written in a form that allows similar renor- 
malization as for the Hamiltonian. Thus, consider the 
identities: 

exp(i q- Ti) = cxp(i q- rij/2) cxp(i q- Rij) (A. 3a) 



exp(iq- fj) = cxp(— ig • rij/2) cxp(i (f • Rij). (A. 3b) 

Using the usual expansion of the exponentials in spherical 
Bessel functions and spherical harmonics, one can show 
that 
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y ^^^Vzf + (-l)')(W,rO|LO)j,(gr,;,/2)jX9i?..)[i1(n.O ® (A.4) 

w 

In the hmit of zero momentum transfer, we indeed find that the electric multipoles can be written as sum of terms 
which factorize the relative and CM contributions of nuclcon pairs. In particular, for L = 2, we regain Eq. HA.2() . 

An important issue in shell model calculations is separation of intrinsic and CM excitations. Although our wave 
functions factorize exactly in intrinsic and CM contributions, one might pick up spurious contributions if one does 
not use translationally invariant operators. In the following, we concentrate on the quadrupole operator. Because 
the last term in Eq. (|A.2|I cannot be summed in a part involving only the A-body CM, one might expect that the 
CM can introduce spurious contributions. We point out, however, that the CM of individual pairs can be changed 
without changing the CM of the A-body system. Furthermore, we can show that there is no contribution from the 
A-body CM. Thus, we start with the translationally invariant expression of the quadrupole operator 

E2 = J20in~RcM), (A.5) 

i 

SO that Eq. (|A.2|l becomes 

0(f, - RcAi) + 0(r, - Rcm) - ^O(ry) + 20(11,, - Rcm). (A.6) 
This can be further transformed by means of Eq. (35) in Ref. [2^ : 
0{n ~ Rcm) + 0{rj - Rcm) = 

\o{nj) + 20{R,j) + 20{Rcm) - ^-^[i?.yYi(i?,,) ® RcmYi{Rcm)]2- (A.7) 

Because = RijYi{Rij) one can sum contributions from all pairs (i, j), so that Eq. HA.7|I becomes: 



0{n - Rcm) = 0{n) + RImY2{Rcm) 

i i 



V47r5! 



A{A- 



l)Rl,,[YiiRcM) 



<E)Yi{RcM)h 



(A., 



Therefore, when computing E2 transition strengths 
with wave functions which factorize the intrinsic and 
Ohfl CM pieces, there are no spurious contributions even 
though the operator usually employed, Eq. IjA.ip . is 



not translationally invariant. The latter differs from its 
translationally invariant form by a term which contains 
only an irrelevant tensor contribution from the CM of the 
A-body system. 



[1] S. Okubo, Prog. Theor. Phys. 12, 603 (1954). 

[2] J. Da Providencia and C. M. Shakin, Ann. of Phys. 30, 

95 (1964); K. Suzuki and S.Y. Lee, Prog. Theor. Phys. 

64, 2091 (1980); K. Suzuki, Prog. Theor. Phys. 68, 246 

(1982); K. Suzuki and R. Okamoto, Prog. Theor. Phys. 

70, 439 (1983). 
[3] K. Suzuki, Prog. Theor. Phys. 68, 1999 (1982); K. Suzuki 

and R. Okamoto, Prog. Theor. Phys. 92, 1045 (1994). 
[4] P. Navratil, J. P. Vary, W. E. Ormand, and B. R. Barrett, 

Phys. Rev. Lett. 87, 172502 (2001). 
[5] P. Navratil, J. P. Vary, and B. R. Barrett, Phys. Rev. 

Lett. 84, 5728 (2000). 
[6] P. Navratil, J. P. Vary, and B. R. Barrett, Phys. Rev. C 

62, 054311 (2000). 



[7] H. Kamada, et. al., Phys. Rev. C 64, 044001 (2001). 
[8] B. R. Barrett, B. Mihaila, S. C. Pieper, and 

R. B. Wiringa, Nucl. Phys. News 13, No. 1, 17 (2003). 
[9] P. Navratil, M. Thoresen, and B. R. Barrett, Phys. Rev. 

C 55, R573 (1997). 
[10] R. B. Wiringa, V. G. J. Stoks and R. SchiaviUa, Phys. 

Rev. C 51, 38 (1995); B. S. Pudliner, V. R. Pandhari- 

pande, J. Carlson, S. C. Pieper and R. B. Wiringa, Phys. 

Rev. C 56 1720, (1997); R. B. Wiringa, Nucl. Phys. A 

631, 70c (1998); S. Pieper and R. B. Wiringa, Annu. 

Rev. Nucl. Part. Sci. 51, 53 (2001). 
[11] R. Machleidt, F. Sammarruca and Y. Song, Phys. Rev. 

C 53, R1483 (1996); R. Machleidt, Phys. Rev. C 63, 

024001 (2001). 



9 



[12] P. Navratil and E. W. Ormand, Phys. Rev. C 68, 034305 
(2003). 

[13] D. C. J. Marsden, P. Navratil, S. A. Coon and B. R. Bar- 
rett, Phys. Rev. C 66, 044007 (2002). 

[14] S. C. Pieper, K. Varga, and R. B. Wiringa, Phys. Rev. 
C 66, 044310 (2002). 

[15] A. C. Hayes, P. Navratil, and J. P. Vary, Phys. Rev. Lett. 
91, 012502 (2003). 

[16] P. Navratil, H. Geyer, and T. T. S. Kuo, Phys. Lett. B 
315, 1 (1993). 

[17] M. Moshinsky, The harmonic oscillator in modem 
physics; from atoms to quarks (Gordon and Breach, New 
York, 1969). 



[18] J. P. Vary, The Many-Fermion Dynamics Code, Iowa 

State University (1992); J. P. Vary and D. C. Zheng, 

ihid. (1994) (unpublished). 
[19] I. Stetcu, B. R. Barrett, P. Navratil, and C. 

W. Johnson, submitted to Int. J. Mod. Phys. E 
arXiv:nucl-th/0409072 . 
[20] F. Ajzcnbcrg-Selove, Nucl. Phys. A490, 1 (1988). 
[21] J. Engel and P. Vogel, Phys. Rev. C 69, 034304 (2004). 
[22] D. A. Varshalovich, A. N. Moskalev, and V. K. Kher- 

sonskii. Quantum Theory of Angular Momentum (World 

Scientific, Singapore, 1988), pp. 167. 



